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ABSTRACT 



We have performed a detailed analysis of the electron temperature and den- 
sity in the the Ring Nebula using the calibrated HST WFC3 images described in 
the preceding paper. The electron temperature (T e ) determined from [N II] and 
[O III] rises slightly and monotonically towards the central star. The observed 
equivalent width (EW) in the central region indicates that T e rises as high as 
13000 K. In contrast, the low EW's in the outer regions are largely due to scat- 
tered diffuse Galactic radiation by dust. The images allowed determination of 
unprecedented small scale variations in T e . These variations indicate that the 
mean square area temperature fluctuations are significantly higher than expected 
from simple photoionization. The power producing these fluctuations occurs at 
scales of less than 3.5 xlO 15 cm. This scale length provides a strong restriction 
on the mechanism causing the large t 2 values observed. 
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Subject headings: Planetary Nebulae: individual (Ring Nebula, NGC 6720) 



1. Background and Introduction 

In the preceding paper we have presented calibrated emission-line images of the Ring 
Nebula obtained with the Hubble Space Telescope's WFC3. These were combined with data 
from ground-based telescope high velocity resolution spectra to develop a detailed 3-D model 
of this prototypical planetary nebula. In this paper we present an analysis of the physical 
conditions. 

The physical conditions in NGC 6720 have been the subject of numerous studies. In 
terms of the spatial coverage, the most complete investigation is that of Lame & Pogge 
(1994), where a Fabry- Perot spectrophotometer was used to map the nebula in the Ha, H/3, 
[N II] 658.3 nm + 575.5 nm, [S II] 671.6 nm + 673.1 nm, [O I] 630.0 nm, [O III] 500.7 nm, 
and He II 468.6 nm lines. Our present study can be viewed as an extension of the Lame 
& Pogge (1994) and Guerrero et al. (1997, henceforth GMC) investigations, but we have 
also observed the 436.3 nm line of [O III] and with more than 20x improvement in spatial 
resolution over the ground-based studies, although with the sacrifice of coverage of about 
one-half of the Main Ring of the nebula. Among the numerous slit spectroscopy studies 
the most recent and comprehensive high resolution investigation was that of O'Dell et al. 
(2007b, henceforth OSH), and the lower resolution multi-line studies of Garnett & Dinerstein 
(2001) and OHS have also appeared. In the preceeding paper we develop and present the 
most up-to-date 3-D model for the nebula. 



In the first part of this paper we characterize the extinction (§ 2.1), the variation in 



the electron temperature (T e ) (§ 2.2[ ) and electron density (n e ) (§ 2.3) in a quadrant of the 



object, derive the Equivalent Widths (EW) of the underlying continuum (§ 2.4), and then 



evaluate the magnitude and importance of small-scale T c variations (§ 2.5) 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under 
NASA Contract No. NAS 5-26555. 

2 Based on observations at the San Pedro Martir Observatory operated by the Universidad Nacional 
Autonoma de Mexico. 
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2. Analysis of the Ring Nebula Calibrated Images 

2.1. Derived Extinction Values 

We have used the calibrated emission-line images (Figure [TJ panels A-E) to determine 
the spatial distribution of the extinction-corrected emission. These images were then used 
to determine the ionization structure within the quadrant of the Main Ring of the nebula 
and the electron temperature and density distributions. 

The extinction correction was determined from the Ha/H/3 line flux ratio. We adopted 
an intrinsic ratio 2.87 based on the predictions for a low density gas at an electron temper- 



ature of 10000 K (Osterbrock & Ferland 2006) and the extinction curve for a ratio of total 



to selective extinction of 3.1 from Table 7.1 of Osterbrock & Ferland (2006). The results are 
best expressed in terms of cnp, the logarithm of the extinction (base 10) at 486.1 nm and 
the results are shown in Figure [l] (panel F). Profiles of cup along the X-axis (PA = 240°, the 
major semi-axis of the elliptical Main Ring) and the Y-axis (PA = 330°, the minor semi-axis) 
are shown in Figure [2j We see there that cup has a radial dependence, being about 0.12 in 
the central core out to the the inner Main Ring, then dropping near the ionization boundary. 



We will see in § 2J2 that there is not a significant variation in T c with distance from the 
Central Star except for a slight increase within 20" of the Central Star. Since 20" is about 
where cjjp begins to drop, it is worth considering if the calculated cup is due to a decrease 
in T e . Osterbrock & Ferland (2006) Table 4.4 shows that the Ha/H/5 ratio decreases by 4 
% in going from 10000 K to 20000 K. The indicated rise in T c towards the center would 
produce a decrease in the Ha/H/5/ ratio of less than 2 %, which would mean that we have 
under-estimated cup by less than 0.025 and the true cup in the core would be even larger 
than is shown in Figure [2j To explain away the outwards decrease in cnp would require T e 
to increase within the Main Ring, but the derived electron temperatures do not show this 
and the outward decrease is likely to be real. 



The average cnp is even higher in regions containing the knots seen in silhouette (O'Dell 
etalpOO^ , which are not numerous in the two samples used in making Figure M The inner 
Main Ring extinction of cjjp = 0.12 is similar to 0.15 found by Garnett & Dinerstein (2001) 
along a line with PA = 91°, and GMC, which found cup = 0.14. The Fabry-Perot study 



(Lame & Pogge 1994) does not give extinction values explicitly, but cite an overall Ha/H/3 
line flux ratio of 3.74±0.22 and 3.67 for the inner 15". These ratios would correspond to 
cnp values of 0.33 and 0.31 respectively. Since these inferred cup values are much higher 
than those from the several other studies, it appears likely that the flux calibration is not 
accurate. 
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Fig. 1. — This mosaic shows the 62.9"x69.4" section of the NGC 6720 images common 
to both the full field and quadrant filters. Panels A through E show the emission-line 
linear images corrected for contamination, if any, for the underlying continuum and non- 
targeted emission-lines, normalized to the maximum features being white. Panel F shows 
the extinction index cnp derived from the surface brightness in Ha and H/5 with the scale 
indicated in Figure [2j The parallel faint linear features passing through the Central Star 
in the cup image are residual "scars" caused by small photometric errors near the edges of 
the two CCD detectors of the WFC3-VIS. The pointing was changed between two otherwise 
identical exposures to allow coverage of the gap between the detectors. The vertical axis is 
pointed towards PA=330°. 
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10 20 30 40 

Distance (arcseconds) 

Fig. 2. — The extinction factor ch/3 is shown as derived from the observed Ha/H/3 ratio and 
assuming an intrinsic ratio of 2.87, which applies to a low density gas at T c = 10000 K. The 
X-axis lies along PA = 240° and the Y axis along PA = 330°. The sample was 4" wide and 
the data have been averaged along bins of ten pixels, corresponding to 0.4". The distance 
shown is from the position of the Central Star. The data closest to that star have been 
edited out because of contamination, thus causing a central gap or in later similar images, 
artificial values near the Central Star. 
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2.2. Electron Temperatures 

We derived the electron temperature for both the [0 III] and [N II] emitting regions using 
equations 5.4 and 5.5 of Osterbrock & Ferland (2006) in their low density approximations. 
The nebular to auroral line ratios were obtained from our emission-line images as described 
in Appendix A. The results are shown in Figure [3] for a 4.0" wide sample along our X axis 
(PA = 240°) and in Figure [4] for a similar sample along our Y axis (PA = 330°). The 
range of distances was limited to those where the relevant surface brightnesses were high. 
[N II] emission arises from only in the Main Ring, giving a nearly constant temperature of 
about 10000 K, with indication of a slight increase at the inner boundary of the Main Ring. 
[O III] is strong from the center of the nebula out to close to the outer boundary of the 
Main Ring. The [O III] temperature is very similar to the [N II] temperature where both 
are available (essentially in the Main Ring). However, the trend suggested by the [N II] 
temperatures is confirmed and extended by the [O III] temperatures. There is a clear rise 
in the [O III] temperature to a central region value of about 11700 K. This inner- region is 
counter-intuitive since photoionization theory indicates that the electron temperature should 
rise with increasing distance from the ionizing star due to the effects of radiation hardening. 
The resolution of this riddle is that either processes affecting the heating become more 
important at the lower densities that apply to the central region or that the central [O III] 
emission arises from material physically further away from the Central Star and is only seen 
in projection at smaller angular separation. 

Electron temperatures have been derived many times for the Ring Nebula, largely in 
connection with efforts to determine elemental abundances. The observational line ratios 
used in their determination have probably improved over the decades and at the same time 
there have been changes in the adopted values of the atomic coefficients such as collision 
strengths and transition probabilities. Limiting ourselves to the most recent work, we note 
that GMC find average values in the Main Ring of 10800 K from [N II] and 11700 K from 
[O III], with central region temperatures of 12600 K from [O III]. Garnett & Dinerstein 
(2001) found the [O III] temperature to be 10000 K in the Main Ring, rising to about 
12000 in the central region. Lame & Pogge (1994) did not detect a variation in their [N II] 
temperature and found a global average of 9643±230 K. Liu et al. (2004) found T ([0 III]) 
= 10600 K and T e ([N II]) = 10200 K. Our temperatures are within the range of these earlier 
spectroscopic results and we confirm the presence of a central increase (which is also seen in 
the Kupferman (1983) study). 
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Fig. 3. — The T c values derived from both the [N II] and [0 III] lines are shown for the 
X-axis profile. The range of results for the two ions was constrained to those where the 
surface brightness of the nebular transitions were the highest. 
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Fig. 4. — Like Figure [3] except for the Y-axis. 
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2.3. Electron Densities 



The densities in the [S II] emitting regions were determined by adopting the relation 
log n e = 4.705 -1.9875 S(671.6)/S(673.1) from Figure 5.8 of Osterbrock & Ferland (2006), 
which is quite accurate over the line ratio range of 1.35 to 0.65, which corresponds to densities 
of about 100 cm~ 3 to 3000 cm~ 3 . The results are shown in Figure [5] for the region of high 
[S II] surface brightness. We see that in the [S II] bright regions of the Main Ring that 
peak values of about 550 cm -3 are found but regions about 350 cm -3 are present. Aller 
et al. (1976) measured [S II] in three points in the Main Ring and when their observed 
line ratios are interpreted with the up-to-date atomic parameters in Osterbrock Sz Ferland 
(2006) yield densities of 600, 360, and 170 cm" 3 . Lame & Pogge (1994) found about 700 
cm -3 for the brightest part of the minor axis (where we find 570 cm _3 )and 500 cm -3 for 
the rest of the nebula. The GMC densities are higher than our own in that they found 615 
cm -3 for their OR sample (radially averaged over 25"-35", which is not the same region we 
have sampled, even though Figure [5] gives results for those distances.). Wang et al. (2004) 
determined densities in the Main Ring from multiple doublets and found 525 cm -3 from the 
[S II] doublet. Given the large local variations in the true densities, we find good agreement 
in the results from the several studies and our derived densities. 



2.4. Equivalent Widths 

The continuum radiation in the Ring Nebula is unusually strong. Since the continuum 
primarily arises from ionized gas processes, it is customary to express the strength relative 
to the the H/3 emission-line, the Equivalent Width (EW). EW is calculated as S(H/3)/Sa, 
where S(H/3) is the H/3 surface brightness in energy units (e.g. ergs) and Sa is the surface 
brightness brightness of the continuum in energy units per wavelength interval (e.g. ergs/A) 
at the same wavelength with the result being expressed in wavelength units (A in our case). 
EW is how wide a sample of the continuum is required to produce as much energy as the 
reference emission-line. This means that increasing strength of the continuum relative to H/3 
produces a smaller EW, an awkward but traditional measure. 

The several atomic processes producing the atomic continuum are detailed in Osterbrock 
& Ferland 2006 and these include for hydrogen free-bound, free-free, and two-photon emission 
and the same processes for contributions from both stages of ionized helium (with the two- 
photon component being relatively unimportant near the wavelength of H/3). The hydrogen 
2S level from which two-photon emission arises can be depopulated by collisions or radiation, 
with the critical density at which the two processes occur equally being 15500 cm" 3 . Since 
we show in § 2.3 that the highest densities are about 500 cm -3 , collisional quenching of 
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Fig. 5. — The electron densities derived from the [S II] 671.6 nm and 673.1 nm lines for an 
assumed electron temperature of 10000 K are shown form both the X-axis (light line) and 
Y-axis (heavy line) profiles. 
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two-photon emission is unimportant. For a pure hydrogen gas at a 10000 K, EW would be 
1400 A. The helium contribution to the continuum is primarily from doubly ionized helium, 
with the emissivity per ion being about seven times that of hydrogen. If the relative number 



abundance of helium to hydrogen is 0.12 and the helium is 70 % doubly ionized (Garnett & 



Dinerstein|2001 ), the expected EW drops to 935 A. For these same conditions the EW drops 
with increasing electron temperature, being 835 A at 12000 K, 770 A at 14000 K, and 750 
A at 16000 K. 

We have determined the EW at H/3 using the properties of the WFC3 filters and adopt- 
ing a constant flux ratio (f A4861 A A5470 = 1-0), finding EW(A) = 690x (Rf487n/Rf547m) -62. 
Figure [6] shows the derived EW along the primary axes of the nebula. If the flux ratio is 
0.90±0.14, as determined from the spectra, then the derived EW values would be about 10% 
larger. The observed WFC3 EW values indicate a central region temperature of about 13000 
K. As an independent check on the unexpectedly low value of EW, we derived the EW from 
the spectra and show the results in Figure [7j Although derived with rigor, there may be a 
small overestimate of the continuum strength (hence smaller EW) due to broad instrumental 
scattering of the very strong [O III] 495.9 nm line. We see good agreement from these two 
methods, with the spectroscopic results being slightly lower, probably because of the instru- 
mental scattered light. In this derivation of EW we noted that there was an almost constant 
background signal in both the F487N and F547M signals, being 2 % of the Main Ring signal 
in F487N and about 10 % for the F547M signal. These backgrounds were subtracted before 
the EW was calculated. 

The surprise is that the EW drops dramatically in the outer parts of the Main Ring. The 
relative continuum is increasing in strength as the outer ionization front is approached. The 
same trend is seen in Figure |8j where in addition to the outer Main Ring drop in EW, there 
are multiple arcs of lower EW values. Although these resemble the arcs of knots in Figure [8j 
the effects of small-scale extinction are nulled-out at the two similar wavelengths used in the 
derivation. The low EW arcs systematically occur within or just outside of the large-scale 
and small-scale [N II] features that indicate He°+H + zones that arise on the Central Star 
side of ionization fronts. This suggests that the lower EW values are due to light scattered 
from the dust particles that must be found in the dense PDR's lying just outside of the 
ionization fronts. This interpretation is strengthened by the fact that the lower EW values 
in the image coincide with the H2 peaks seen in Figure 3 of the preceding paper. 

In order to quantify the amount of excess continuum emission, we constructed templates 
from the H/3 and He II line images to model the atomic continuum from H + and He ++ , 
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Fig. 6. — The Equivalent Width of the H/3 emission line and its underlying continuum 
are shown along the PA values 60°, 150°, 240°, and 330° as determined from the WFC3 
F547M and F487N images. Decreasing Equivalent Width values indicate relatively stronger 
continuum values. Most of the regions contaminated by starlight have been edited out. 



13 




Fig. 7. — Like Figure [6j except that the EW values are derived from the ground-based 
spectra used for calibration. Different PA values are indicated. The regions contaminated 
by starlight have been edited out. 
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Fig. 8. — This image of the derived Equivalent Width of the H/3 emission line is for the same 
FOV as Figure 3 of the preceding paper. The stitching caused by the dithering across the 
gap between the two CCD's is even more obvious, but does not affect the derived equivalent 
width. The positive Y axis points towards PA = 335°. The variations in the Equivalent 
Width are discussed in the text. 
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respectively]^] The gas temperature (assumed constant within the nebula) was adjusted in 
order that the predicted atomic continuum be equal to the observed continuum in regions 
of the nebula with the largest EW). The resultant value was 14,000 K, which is significantly 
higher than the temperatures determined from the [O III] and [N II] lines. Photoionization 
models confirm that the highly ionized He ++ +H + region should be hotter than regions where 
helium is singly ionized ( Osterbrock & Ferland 2006 ) and the high kinetic temperature needed 
to account for the H/3 equivalent width seems reasonable. The modeled atomic continuum 
is then subtracted from the observed continuum to give the residual continuum outside of 
the inner region where EW is a maximum. This residual emission is found to be very well 
correlated with the H2 2.12 /im emission (Figure 3, Panel F of the preceding paper), which 
supports the view that it is produced in neutral/molecular gas and is due to dust scattering. 

The total flux of residual nebular emission at 5470 A calculated by this method (after 
removing all stars in the image and correcting for an extinction of cnp = 0.15) is 2.45 x 
10~ 14 erg s _1 cm -2 A -1 , which is roughly 25% of the atomic continuum from the nebula. 
For comparison, the total flux from the nebula's central star is 2.06xl0~ 15 erg s _1 cm -2 
A -1 , which is over 10 times less. This means that the central star cannot be the source of 
the scattered continuum seen in the nebula, and the illumination must come from outside. 



The average intensity of the interstellar radiation field at this wavelength (Black 1987) is 
1.3xl0 -7 erg s -1 cm -2 sr _1 A -1 and, since the Ring Nebula is about 170 pc above the 
Galactic plane, the illumination is from one hemisphere only, so that the angle-averaged 
mean intensity, J\ is half this value. If the scattering optical depth of the dust in the nebula 
is r < 1, then the required value of J\ in order to produce a scattered flux F\ at the earth 
is J\ ~ Fx/ (rfl), where fl is the solid angle subtended by the nebula. Taking an upper limit 
of r ~ 1 and an average radius of 45" for the nebula requires an external radiation field 
of Jj, >1.6xl0~ 7 erg s _1 cm -2 sr _1 A -1 , which is roughly two times larger than expected 
for the Galactic diffuse field. The agreement is fair considering the uncertainties involved, 



especially since the Galactic field is expected to vary from place to place (Parravano et al. 
2U031. 



3 Since we have no WFC3 image of an He I line, we also use H/3 as an approximate proxy for the He + 
continuum. This is justified since the He + continuum is a small fraction (roughly 5%) of the H + continuum 
and the ground-based spectra show that He I/H/3 varies very little within the nebula (see Fig. 16). 
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2.5. Small Scale Temperature Variations, the t 2 Problem and the Abundance 

Discrepancy Factor 

There is an elusive problem in gaseous nebulae that has to do with their physical condi- 
tions. The temperatures derived from collisionally excited lines (CL) are higher than those 
derived from recombination lines (RL) and from Balmer continua. The difference is consid- 
erably higher than that predicted by sophisticated photoionization models. Moreover it has 
important implications for the derivation of chemical abundance ratios. The abundances 
derived from collisionally excited lines assuming thermal homogeneity are typically about a 
factor of two to three times smaller than those derived from recombination lines. To reconcile 
both types of abundances it is necessary to invoke the presence of temperature variations 
considerably higher than those predicted by photoionization models. 

We address the problem of finding the mechanism or mechanisms responsible for these 
temperature variations by studying the distribution of the temperature across the face of the 
Ring Nebula based on the very high spatial resolution provided by the HST observations. 
The filters available permit us to determine highly accurate T(0 ++ ) and T(N + ) values. 
There are two observational quantities that are most useful: the determination of the mean 
square temperature variation in the plane of the sky and the study of the typical spatial size 
of the temperature variations. In order to obtain the observed mean square temperature 
variation in the plane of the sky we will follow the set of equations introduced by Peimbert 
(1967) and O'Dell et al. (2003). 

The most meaningful measure of small scale temperature variations is the volumetric 
variation, for which we define the pair of equations for the average temperature, T (X +t ), 
and the mean square temperature fluctuation, t 2 (X +t ), given by 

T(v+ , ] _ jT e (r)n e (r)n(X^;r)dV 
01 ' Jn e (r)n(X+ i ;r)dV ' {) 

and 



,2, Y+i , _ J(T e -T (X^n e n(X +i )dV 
1 j " Tq j n e n(X +i )dV ' [) 

respectively, where T e and n e are the local electron temperature and density, n(X +l ) 
is the local ion density corresponding to the observed emission line, and V is the observed 
volume. 



However, our observational results are limited to determination of variations in columns 



-17- 



across the plane of the sky t\(0 ++ ) and t\(N + ). This requires the definition of a columnar 
temperature T c 

T (X^ a 6)- f^MX+^dl 

Tc{X fn e n(X +i )dl ' ( ) 

where a is the right ascension and 5 is the declination of a given line of sight (corre- 
sponding to a given pixel). Then, the average temperature can be written as 

rp * y+i , _ jT c (X^a,5)Jn e n(X +l )dldA 
01 ] J J n e n(X +i )dldA ' K ) 

and the mean square temperature variation over the plane of the sky, t\(X +l ), can be 
defined as: 

,2 ( Y+i) _ J(Tc(X«) - UX^)f J n e n(X +i )dldA 
A{ ] T {X+ i ) 2 J J n^X+^dldA ' K ' 



2.5.1. Determination oft 2 A (0 ++ ), and t\(N + ) from the WFC3 images 

In order to isolate the effects of the small-scale temperature fluctuations in the plane 
of the sky, we first "detrend" the temperatures by removing the systematic radial variation. 
For [O III] we adopted a mean temperature that varies with projected radius, R, from the 
central star as 

T (O+ + , R) = ^(T in + T out ) + ^(T out - T in ) tanh (^j^j , (6) 

with parameters determined from fitting to the observed temperatures: inner temperature 
T in = 11,350 K; outer temperature T out = 9560 K; transition radius R = 17.6"; transition 
width 5R = 4.4". In the case of [N II], no clear radial trend is visible (Figure [3] and 
Figure [4]) so a constant mean temperature was used. The uncertainties in the pixel-to-pixel 
temperature measurements were dominated by uncertainty in the measurement of the auroral 
lines, [N II] 575.5 nm and [O III] 436.3 nm, since these were roughly 100 times fainter than 
the corresponding nebular lines. The principal contributions to these uncertainties were 
(1) shot noise associated with the Poisson statistics of the arriving photons; (2) readout 
plus dark noise associated with the CCD; (3) pixel-to-pixel flat-field variations; (4) cosmic 
rays artifacts. Noise sources (2) and (3) are very well understood and characterized for the 
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WFC3 instrument, and their effect on the measured t 2 was found to be negligible. Source 
(1) is also well understood and will produce a contribution to the observed fluctuations of 
^Lise = 1 / ' i. 3 " 2 'N) j where N is the number of detected photoelectrons and s = d\iaT/d\nr is 
the sensitivity of the derived temperatures to the diagnostic line ratio r (at 10 4 K s = 2.5 
for [N II] and s = 3.29 for [O III]). Source (4) was the most problematic, since 5-9% of the 



pixels in each WFC3 exposure are affected by cosmic ray artifacts (Rajan et al. 2010), and 



considerable care had to be taken taken to minimize its effects on our analysis. 

Multiple exposures were taken of the auroral lines (3 exposures of 900 s each for the 
FQ436N filter, 4 exposures of 600 s each for the FQ575N filter), which were combined in the 
standard pipeline reduction process that carries out automated cosmic ray artifact rejection 



(Bushouse 2010). Although the automated pipeline process effectively removes the bright 
cores of cosmic ray artifacts, many cosmic ray artifacts possess fainter halos that are still 
visible on the combined images. We therefore employ two filters in order to restrict our t 2 
analysis to only the highest quality portions of the image. First, we rejected all pixels which 
were flagged as bad by the pipeline process in even one exposure. Second, we rejected all 
pixels in which the standard deviation between the three (or four) individual exposures was 
larger than a certain factor times the expected width of the Poisson distribution, yN. We 
found that a value of 1.5 for this factor is a good compromise between effectively removing 
low-level cosmic ray artifacts, while at the same time not throwing away too many pixels. 
The fraction of image pixels that are rejected by these two filters was 30 to 40%. 

In order to both reduce the effects of Poisson noise and investigate the spatial scale of 
the plane-of-sky temperature fluctuations, we calculated t\ for a series ofmxm binnings of 
the images: lxl pixel, 2x2 pixels, 4x4 pixels, etc., up to 128 x 128 pixels. The binning 
was performed by averaging the surface brightnesses of the auroral and nebular lines, but 
for only the good pixels (defined as in previous paragraph) in each m x m patch. The mean 
temperature of each patch was calculated as in § 3.3 and t\ was calculated as in Equation (7), 
but using the surface brightness of the nebular line as a proxy for J n e n(X. +i ) dl. At the 
same time, the sum of the total number of detected photoelectrons that contribute to each 
patch was tracked, so that t 2 oise could be calculated as above. 

The results are shown in Figure [9] for [O III] and [N II] respectively. The cross symbols 
show the measured t\, while the solid line shows the expected contribution from noise, 
which, in addition to Poisson and read-out noise, also includes the a very small (< 10~ 4 ) 
contribution from flat-field variations. The plus symbols are the result of subtracting the 
expected noise contribution from the observed values. It should be noted that the diffraction- 
limit resolution of the images is about 2 pixels and so the results of binning at 1 x 1 pixels 
should be considered lower limits on the true fluctuations at that scale. 
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Fig. 9. — Global plane-of-sky temperature fluctuations t\ for (a) [O III] and (b) [N II] as a 
function of bin area (lower horizontal axis) or angular scale (upper horizontal axis). Cross 
symbols show the measured values, while plus symbols show the values after correcting for 
the expected contribution from noise (indicated by the solid line). The left-most point in 
each graph is affected by the finite angular resolution of the telescope and hence is only a 
lower limit to the true t\ at that scale. In the case of [O III], "detrended" temperatures are 
used, from which the large-scale radial variation has been removed as described in the text. 
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The noise-corrected results for both ions can be well fit by piecewise power laws of the 
form t\ oc L a , where L is the angular scale, with a break in the power index a at a scale of 
-^break — 0.5" as shown in Table 1. The fall off of the fluctuations with scale is significantly 
steeper at scales smaller than Lbreak ( a — — 5 f° r [O HI] and a ~ — \ for [N II]) than at larger 
scales (a ~ — | for [O III] and a ~ — ^ for [N II]). 



2.5.2. Determination of the t 2 volumetric factor for selected ions 

Up to this point we have presented the results in terms of t\(0 ++ ) and t^(N + ), but 
to compare with photoionization models and with observations of the whole object what 
is needed are the variations in three dimensions t 2 {X +l ). In order to obtain the total 
t 2 (X + *)values we need to consider the variations along the line of sight. It can be shown 
that the relevant equation is: 

t 2 (X +i )=t 2 A (X^) + (t 2 c (X^,a,5)), (7) 
where t 2 , the variation along a given line of sight, is given by: 



t 2 (X^ a 5)- I^n(X^)(T e -T c ) 2 dl 



and the average over all lines of sight is given by: 



, „ xN f t 2 T 2 f n e n(X +i )dldA 

Liu et al. (2004a,b) have presented slit spectra line intensity observations of the Ring 
Nebula that give sample regions across the whole object, while our HST observations include 
approximately one quadrant. We consider our highly sampled quadrant to be representative 
of the whole object. From the data of Liu et al. it is found that O is distributed as follows: 
32.4% as + , 52.3% as ++ and 15.3% of O in higher degrees of ionization. The last number 
comes from the assumption that the ratio of abundance of the higher degrees of ionization 
to n(0 + + ++ ) is given by n(He ++ )/n(He + + Hc ++ ) and that the He/H abundance ratio 
is 0.112. 

We can make three estimates of Tq(0 ++ ) and t 2 (0 ++ ) from the data by Liu et al.(2004a,b). 
a) By assuming that the nebula is chemically homogeneous and that the O II RL and the 
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[O III] CL originate in the same volume we can obtain two temperatures for the ++ re- 
gion, one given by the [O III] 4363/4959 ratio and the other by the O III 4959 to the O II 
VI multiplet ratio (these two ratios have very different temperature dependences) and by 



using the formalism presented by (Peimbert 1967), and the computations by Storey (1994) 



( Bastin fc Storey||2006 ) and Storey (see Liu (2012)) we have determined the T and t 2 values 
presented in Table 2. b) By assuming that the average temperature and the temperature 
variations are similar in the + and ++ regions we can combine T(4363/4959) with T(He 
1,5876/7281) and obtain the To and t 2 values presented in Table 2. c) By assuming that 
T(4363/4959) is representative of the whole object we can compare it with T(H11/Balmer 
Jump) and obtain the To = and t 2 values in Table 2. The three t 2 values presented in Table 
2 are in good agreement and a representative t 2 (0 ++ ) value is in the 0.045 to 0.055 range. 

It is not possible to derive t 2 and To combining T(Hel) and T(H11/Balmer Jump) 
because the temperature dependence of both determinations is almost the same. 



An absolute lower limit of the t\(0 ++ ) value derived in 2.5.1 is 0.009, this result is 
obtained by sampling all the observed area with bins of 2x2 pixels. The true value of 
t\(0 ++ ) could only be obtained with infmitesimally small sampling. There is a clear trend 
of increasing t\(0 ++ ) values with smaller linear scale, suggesting that the true value is 
higher; but clearly this trend cannot continue to arbitrarily large values. By extrapolating 
the t 2 A (0 ++ ) values for 8x8, 4x4 and 2x2 pixels to a sampling area of lxl pixel a t^(0 ++ ) 
value of 0.013 is obtained (see Figure |9|. This value should be considered as a soft lower 
limit to the true t\(0 ++ ) value since the value could be higher if this trend continues to 
a sampling with areas smaller than lxl pixel. The 0.013 value is about a factor of four 
smaller than the volumetric t 2 (0 ++ ) value (see Table 2). This difference is expected, since 
the temperatures for each point in the plane of the sky are an average of the temperatures, 
see equation (5) and equation (10). As such t\(0 ++ ) is only a lower limit to t 2 (0 ++ ). The 
t\/t 2 ratio is also expected to be small when the characteristic scale of the spatial variations 
is much smaller than the whole object. 

From geometrical considerations one would expect the variations of one additional di- 
mension, t 2 , to be at least half as large as the variations of two dimensions, t 2 A \ this works 
for large scale temperature variations. The power included in t 2 increases in the presence of 
small size scale variations. This is because the many independent thermal elements along 
each line of sight will be masked in the averaging process, greatly reducing the total t\ 
compared to t 2 . 

Similarly, an absolute lower limit of t\(N + ) value derived in the previous subsubsection 
is 0.0053. This result is obtained by sampling all the observed area with bins of 2x2 pixels. In 
this case the trend, showing larger inhomogeneities at small scales, is also visible. Compared 
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with the t 2 A (0 ++ ) determinations the t^(N + ) shows more power at scales larger than 8x8 
pixels and a slightly shallower behaviour at smaller scales. As is the case for ++ the "most 
important scale" is probably smaller than our detection threshold. By extrapolating the 
t\(N + ) values for 8x8, 4x4 and 2x2 pixels to a sampling area of lxl pixel, a t\(N + ) value 
of 0.008 is obtained (see Figure [9]). 

The previous discussion implies that the temperature variations are dominated by varia- 
tions in small distance scales of the order of a thousandth of a parsec or less. At the adopted 
distance of the Ring Nebula the length of two pixels amounts to about 9xl0 14 cm. 

For comparison O'Dell et al.(2003) obtained for the Orion Nebula t 2 A (0 ++ ) = 0.0079 
and t 2 (0 ++ ) = 0.021 ±0.005, a ratio of t 2 (0++) over t 2 A (0 ++ ) of about 2.0-3.3. 

2.5.3. The O/H abundance discrepancy factor 

New Location as first paragraph of this section. The recombination line abundances of 
H, He, C, N, O, and Ne have a similar dependence on the electron temperature, therefore 
in the case of chemical homogeneity the abundances derived from the ratio of two RL are 
almost independent of the temperature structure. This is not the case for the abundances of 
collisionally excited lines of C, N, O, and Ne relative to H because the first four do depend 
strongly on the temperature structure. Their intensity increases with temperature; while the 
opposite occurs with the recombination lines. In addition, their intensity decreases mildly 
with increasing temperature. Consequently in the presence of temperature variations the 
preferred abundances are those derived from the ratio of two recombination lines. However, 
an accurate derivation requires measurement of inherently weak and some-times blended 
lines. This demands high signal to noise spectra of excellent wavelength resolution. 

In Table 3 we present the chemical abundances derived by Liu et al. (2004a,b) for the 
Ring Nebula based on RL and on CL lines. The total O abundances are smaller by 0.04 dex 
than those presented by Liu et al. (2004b) because they did not correct properly for the 
fraction of +3 based on the He + and He ++ abundances. From Table 3 it can be seen that 
the recombination abundances for C, N, O, and Ne are higher than the collisional abundances 
by about 0.4 dex. 

There are in the literature two families of results to explain the CL versus the RL abun- 
dance differences based on the temperature dependence of the abundance determinations: 
a) the presence of low temperature H-poor inclusions that are rich in helium and heavy 
elements in planetary nebulae (e. g. Liu 2006 and references therein), or b) the presence 
of temperature variations in a homogeneous medium due to additional physical processes 
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not included in photoionization models (e. g. Peimbert and Peimbert 2006, and references 
therein). The fact that, T(H), T(He) and T(Vl/4959) are similar indicates that an impor- 
tant chemical inhomogeneity is not present in the Ring Nebula. We are left with the idea 
of a well mixed material with temperature fluctuations. The details of the origin of these 
fluctuations remain unexplained. Physical mechanisms known to play at least a small part 
in explaining these fluctuations are: shocks, magnetic reconnection and shadowed regions, 
the relative importance of these processes must be studied further. 



2.5. 4- Photoionization model predictions 

We computed a simple model of the main H II region of the Ring Nebula using reasonable 
parameters. These included a central star temperature of 1.2xl0 5 K, a luminosity of 200 
L Q , and abundances typical of the Ring Nebula as determined from collisionally excited 
lines. This model reproduced the overall ionization structure and the intensities of the 
stronger lines. We then determined the temperature variations across the model, using the 
methods described in Kingdon & Ferland (1995). The "structural t 2 ", the mean variation 
in temperature across various emission zones, was 0.0062 for the H II region, and 0.0026 
across the [O III]-forming region. This model also predicts for the He ++ emitting zone that 
T c = 12500 K, essentially in agreement with the value of 13000 K inferred from the EW in 
the central region ( ^2.4 ). 



From observations (§ |2.5.2| we have determined that t 2 (0 ++ ) = 0.050±0.008, and t 2 4 (O ++ )>0.0090. 
These values are considerably higher than the t 2 (0 ++ ) = 0.0026 value predicted by the pho- 
toionization model. These two results imply that in addition to photoionization there must 
be one or more additional physical processes responsible for the observed t 2 values. Unfortu- 
nately, although the relevant scales for these processes (< 0.5") are in principle accessible to 
HST observations, the low signal-to-noise of the auroral line images are a serious impediment 
to determining their origin. 



2.5.5. Possible origins of small-scale temperature fluctuations 

In efforts to explain the existence of small-scale temperature fluctuations in photoion- 
ized nebulae, various physical mechanisms have been proposed in the literature. These 
mechanisms can be divided into two classes: first, those that posit an extra source of energy 
in addition to the photoionization heating that is generally taken to dominate the heating 
within the nebula, and, second, those that attempt to modify the heating/cooling balance 
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through other means. 



An example of the the first class of mechanisms is the action of highly supersonic stellar 



winds, either directly through the thermalization of the wind kinetic energy (Luridiana et 



al. 2001), or via conduction fronts at the boundaries of hot, shocked bubbles (Maciejewski 



et al. 1996). In the case of H II regions, quantitative studies have suggested that these 



mechanisms are unable to satisfy the derived energy requirements (Binette & Luridiana 



2000), but we are unaware of comparable calculations tailored to the case of Planetary 



Nebulae. A more promising source of energy is perhaps internal turbulent motions driven by 
transonic photoevaporation flows from dense condensations at the periphery of the nebula 
(such as the knots observed in the Ring Nebula, see § 2.2.1 of the preceding paper, and 
where the dissipation of the turbulent energy would be principally via low Mach number 
shocks. In the case of H II regions, simulations indicate that the total kinetic energy of 



such turbulent motions is similar to the thermal energy of the ionized gas (Arthur et al. 



2011), and similar processes are expected to occur in Planetary Nebulae ( Garcia- Segura et 



al. 2006). Another potential source of extra heating is magnetic reconnection Lazarian & 



Vishniac (1999); Lazarian et al. (2012), although no quantitative estimate of its importance 



in the context of Planetary Nebulae has yet been made. 

An example of the second class of mechanism is the possible existence of small-scale 



pockets of high-metallicity gas (Kingdon & Ferland 1998 Robertson- Tessi & Garnett 2005 



Liu 2006), which, due to vastly enhanced collisional metal line cooling, will have a lower 



equilibrium temperature than the rest of the nebula. The origin of these pockets may be 
solid bodies such as comets, which were evaporated during the asymptotic giant branch 



phase that preceded the formation of the planetary nebula (Henney & Stasihska 2010). 



Another example of this class of mechanism is the existence of diffusely illuminated regions 



of the nebula that are shadowed from the direct stellar radiation (O'Dell et al. 2003). These 



too would be lower temperature than the general nebula, due to being illuminated by the 
relatively soft ionizing field characteristic of recombination radiation. 



3. Summary and Conclusions 

We have used the flux-calibrated HST WFC3 emission-line images described in the 
preceding paper to determine the physical conditions in the Ring Nebula down to sub-arcsec 
scales. Our major conclusions are listed below. 

1. The equivalent width of the continuum using the H/3 emission-line as a reference 
is unexpectedly small, indicating a strong relative continuum. In the inner regions of the 
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nebula this is probably due to the electron temperature being as high as 13000 K and in the 
outer parts of the nebula it probably indicates scattering of starlight in the diffuse Galactic 
radiation field, although a value twice the global average of the diffuse Galactic radiation 
field is required for this interpretation to apply. 

2. The electron temperature as determined from nebula to auroral transitions of singly 
ionized nitrogen [N II] and doubly ionized oxygen [O III] varies systematically. T c increases 
from at or slightly lower than 10000 K in the outer region (from [N II]) to about 11500 K 
in the inner- most region (from [O III]). Near the central star there is a region of He ++ +H + , 
where oxygen will be triply ionized. This means that lines of sight near the central star 
measure material in the Lobes and the higher temperature derived from the equivalent 
width is characteristic of the innermost region of the nebula. 

3. About half of the power in the observed t^(0 ++ ) and t^(N + ) values originates in 
structures smaller than 8 pixels (3.5 x 10 15 cm). After correcting for the presence of structures 
smaller than 2x2 pixels as well as the contribution of the variation of the temperature 
along the line of sight, it follows that most of the volumetric t 2 comes from structures 
with typical sizes of 2x2 pixels (9 x 10 14 cm). These small scales suggest the presence of 
one or more physical processes, small scale chemical inhomogeneities, radiation shadowing, 
magnetic reconnection, or dissipation of turbulent energy in shocks as possible causes for the 
large t 2 values. 
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Table 1. Contributions to noise-corrected plane-of-sky temperature fluctuations 



Source 


t\(\0 III]) 


ti([N II]) 


Large-scale radial gradient (L > 5") 
Fluctuations on medium scales (L = 0.5" to 5") 
Fluctuations on small scales (L = 0.08" to 0.5") 


0.004 
0.001 (oc L-VS)) 
0.004 (oc L~ 4 / 5 )) 


~ 

0.002 (oc L- 1 / 10 ) 
0.003 (oc L- 1 / 2 ) 


Total nebular plane-of-sky fluctuations (L > 0.08") 


0.009 


0.005 



Table 2. T and t 2 Determined from Slit Spectroscopy* 



Temperatures Used 



8775±270 0.057±0.008 T(4363/4959) = 10630 and T(Vl/4959) = 8810 

9050±370 0.052±0.011 T(4363/4959) = 10630 and T(Hc I, 5876/7281) = 8290 
9630±750 0.032±0.022 T(4363/4959) = 10630 and T(H11/Balmcr Jump) = 9100 



Note. — *The line and continuum fluxes are from the entire nebula as described in the 
text. 



Table 3. The Ring Nebula Abundances Determined from Slit Spectroscopy 





H 


He 


C 


N 


O 


Ne 


RL 


12.00 


11.05 


9.10 


8.62 


9.14 


8.71 


CL 






8.59 


8.22 


8.76 


8.23 



Note. — All abundances are in base 10 logarithms 
relative to the abundance of H set to 12.00 



